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A new approach to the production and detection of bound neutron clusters is presented. The technique 
is based on the breakup of beams of very neutron-rich nuclei and the subsequent detection of the recoiling 
proton in a liquid scintillator. The method has been tested in the breakup of 11 Li, 14 Be and 15 B beams by 
a C target. Some 6 events were observed that exhibit the characteristics of a multineutron cluster liberated 
in the breakup of 14 Be, most probably in the channel 10 Be+ 4 n. The various backgrounds that may mimic 
such a signal are discussed in detail. 



1 Introduction 

The very light nuclei have long played a fundamental role in 
testing nuclear models and the underlying nucleon-nucleon 
interaction. Whilst much effort has been expended in at- 
tempting to model stable systems, a number of ambiguities 
remain. For example, the ground states of 3 H and 3 ' 4 He 
do not appear to be particularly sensitive to the form of 
the interaction (!]]. In this context the study of systems 
exhibiting very asymmetric N/Z ratios may provide new 
perspectives on the N-N interaction and few-body forces. 
In the case of the light, two-neutron halo nuclei such as 
6 He, insight is already being gained into the effects of the 
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three-body force |2[]. Very recently evidence has been pre- 
sented that the ground state of 5 H exists as a relatively 
narrow, low-lying resonance ||. In the case of the lightest 
N = 4 isotone, 4 n, nothing is known |Q, The discovery of 
such neutral systems as bound states would have far reach- 
ing implications for many facets of nuclear physics. In the 
present paper, the production and detection of free neutron 
clusters is discussed. 

The question whether neutral nuclei may exist has a long 
and checkered history that may be traced back to the early 
1960s H • Forty years later the only clear evidence in this re- 
spect is that the dineutron is particle unstable. Although 3 n 
is the simplest multineutron candidate, the effects of pairing 
observed on the neutron drip-line suggest that 4 > 6 ' 8 n could 
exhibit bound states Q. Concerning the tetraneutron, an 
upper limit on the binding energy of 3.1 MeV is provided 
by the particle stability of 8 He, which does not decay into 
a+ 4 n. Furthermore, if 4 n was bound by more than 1 MeV, 
a+ 4 n would be the first particle threshold in 8 He. As the 
breakup of 8 He is dominated by the 6 He channel 0, the 
tetraneutron, if bound, should be so by less than 1 MeV. 

The majority of the calculations performed to date sug- 
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gest that multineutron systems are unbound Q]. Interest- 
ingly, it was also found that subtle changes in the N-N 
potentials that do not affect the phase shift analyses may 
generate bound neutron clusters Q. In addition to the 
complexity of such ab initio calculations, which include the 
uncertainties in many-body forces, the n-n interaction is 
the most poorly known N-N interaction, as demonstrated 
by the controversy regarding the determination of the a nn 
scattering length || . The lack of predictive power of cal- 
culations of few-body systems at the 1 MeV level (see for 
example Ref. ||) therefore does not exclude the possible 
existence of a very weakly bound n. 

Experimentally, no limit has yet been placed on the bind- 
ing energy of 4 n (or any other ^n). Rather only limits on 
the production cross-sections could be estimated from two- 
step P, [l| or direct reactions @ [D| |l§ [j/j], 
|l9| , |20| . In this paper, a new technique to produce and 
detect neutron clusters is presented. The method is based 
on the breakup of an energetic beam of a very neutron-rich 
nucleus and the subsequent detection of the ^n cluster in 
a liquid scintillator. This work represents the continuation 
of an experimental programme investigating the structure 
and, in particular, the correlations within two-neutron halo 
systems HJ || |f & 

The paper is organised as follows. Section 2 describes 
the experimental technique and analysis procedures. The 
results, including the observation of some 6 events exhibit- 
ing characteristics consistent with a multineutron cluster 
liberated in the breakup of 14 Be, are presented in section 3. 
A detailed discussion of the results is given in section 4. 
Special attention is paid to the various backgrounds (most 
notably pile-up) that may mimic the signal arising from a 
multineutron cluster. Finally a summary and outlook on 
future work is briefly presented in section 5. 

2 Experimental technique 
2.1 Previous experiments 

There have been essentially two categories of experiments 
that have searched for A n systems. The first consists of 
the production of A n in reactions such as neutron-induced 
fission of U flG, [Ll| or proton and light-ion fragmentation 



of a heavy target |I2| , |13[ . Any recoiling A n are then, in 
principle, signaled by the radiochemical separation of decay 
products from ( A n,xn) reactions in a secondary target. An 
extremely pure target and a detailed analysis of all possible 
backgrounds are thus needed. As such only upper limits 
for the ^n production cross-section, assuming cross-sections 
for the ( A n,xn) reactions, could be determined. The only 
positive claim []l2f was later explained as arising from an 
underestimation of the production of very energetic tritons 
& 

The second class of experiments involves direct reac- 
tions of the type a(b, c) A n, where discrete values of the 
energy of the ejectile, c, correspond to states in the ^n sys- 
tem. This technique can thus also reveal unbound states. 



These searches have relied on the very low cross-section 
(typically ~ 1 nb) double-pion charge exchange (D7rCX) 
[14, [la, |l6| and heavy-ion multinucleon-transfer reactions 
[17, [l^, |l9[ |2(| . Again no conclusive evidence for a bound 
^i or resonant states has been found in these studies. 
There are many problems inherent in this technique: the 
precise knowledge of the many-body (A+l) phase space; 
the background from target impurities; and the bias intro- 
duced by the fact that both the A n and the ejectile have to 
be formed in the reaction, lowering further the production 
cross-section. 

2.2 Principle 

Clustering appears in many light nuclei close to particle 
emission thresholds [^5|. Examples include a+t cluster- 
ing in 7 Li (threshold at 2.5 MeV), a+a+n in 9 Bc (S n = 
1.6 MeV), and a+4n+a in 12 Bc (threshold at 12.1 MeV). 
In light neutron-rich nuclei, components of the wave func- 
tion in which the neutrons present a cluster-like configu- 
ration may be expected to appear pq . Owing to pairing 
and the confining effects of any underlying a-clustering, the 
most promising candidates may be the drip-line isotopes of 
Helium and Beryllium, 8 He (S in = 3.1 MeV) and 14 Bc 
{S 4n = 5 MeV). 

We have investigated existing data for the breakup of a 
35 MeV/N 14 Be beam by a C target ||, ||, Other com- 
ponents present in the beam which will be exploited here 
were n Li at 30 MeV/N and 15 B at 48 MeV/N. In such re- 
actions relatively high cross-sections (typically ~ 100 mb) 
are encountered. Consequently, even only a small compo- 
nent of the wave function corresponding to a multineutron 
cluster could result in a measurable yield with a moderate 
secondary beam intensity. Furthermore, the backgrounds 
encountered in D-7rCX and heavy-ion transfer reactions are 
obviated in direct breakup. The main difficulty in the ap- 
proach lies in the direct detection of the A n cluster. 

The details of the experimental setup have been described 
elsewhere |2^, p8[ . Therefore only the salient fea- 

tures are recalled. Following the breakup, the charged frag- 
ments were detected using position-sensitive Si-Csl tele- 
scopes and the neutrons using 90 modules of the DEMON 
array, located at distances of 3.5 to 6.5 m downstream of 
the target ^3|. The energy of the neutrons (E n ) was de- 
rived from the time-of- flight, the start being furnished by a 
thin Si detector placed just forward of the target. An av- 
erage resolution of 1.5 ns was obtained. For particles with 
A > 2, E n corresponds to the energy per nucleoli. Stan- 
dard pulse-shape discrimination methods based on the light 
output from the liquid scintillator were employed to iden- 
tify the neutrons from the 7 and cosmic-ray backgrounds 
p9|| . Except where noted, only events in which precisely 
one DEMON module fired in coincidence with a charged 
fragment were considered in the analysis, in order to avoid 
any possible contributions from cross-talk |2£| . 

The predominant mechanism for the detection of neu- 
trons in a liquid scintillator such as that used in DEMON 
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FIG. 1: Light output calibration for a DEMON module. The 
open symbols correspond to the Compton edge of 7-rays from 
22 Na (511 and 1275 keV), 137 Cs (662 keV), 60 Co (1333 keV) 
and 241 Am 9 Be (4.44 MeV) sources, and the closed symbol to 
cosmic-ray muons Ji?], ps[ |. The solid line is a parabolic fit. 



is n-p scattering |3(|, in which the proton recoils with an 
energy (E p ) up to that of the incident neutron. In general, 
the neutron does not lose all its energy in the interaction 
and may escape from the detector |28| . The energy of the 
recoiling proton can be determined from careful source and 
cosmic-ray calibrations of the charge deposited in the mod- 
ule |27], |2{| Q (Fig. [I]) . This may then be compared to the 
energy per nucleon of the incident particle derived from the 
time-of-flight (E n ). For a single neutron and an ideal de- 
tector E p /E n < 1. For a real detector the finite resolutions 
can give a higher limit, and for DEMON this is ~ 1.4. In 
the case of a multineutron cluster, E p can exceed the inci- 
dent energy per nucleon and E p /E n may take on a range 
of values extending beyond 1.4, as shown in Fig. |^ — the 
scale on the upper axis indicates the maximum value as a 
function of the multineutron mass number. 



Reactions on Carbon in the liquid scintillator, as for ex- 
ample n-C scattering or C(n,3a), do not present any dif- 
ficulties to the present technique as the associated light 
outputs translate to very low E p /E n p0| . In addition, ex- 
cept for n-C scattering, the cross-sections for reactions on 
C are negligible compared to scattering on Hydrogen in the 
energy range considered here. 



2.3 Calibrations and energy range 

As is evident in the source and cosmic-ray calibrations in 
Fig. [I], the DEMON modules exhibit saturation effects at 
very high light output, as observed in earlier work f29| . It 
should be noted that the initial goal of the experiment run 
to acquire the data analysed here was not the search for 
multineutron clusters, and as such the analysis of very high 




FIG. 2: Distribution of the ratio of proton energy, E p (MeV), to 
the energy derived from the flight time, E n (MeV/N), for data 
from the reaction C( 14 Be, 12 Be+n) (histogram) and for simu- 
lations of elastic scattering of 1,3,4 n (solid, dashed and dotted 
lines, respectively) on protons. The experimental resolution has 
been included in the simulations. 



light outputs was not foresecnf]. The total charge versus 
light output is well described up to ~ 25 MeVee using a 
parabolic adjustment (Fig. [l]), as found in previous tests 
which included measurements with a 15.5 MeV 7-ray p9[ . 
A deposited charge of ~ 25 MeVee corresponds to a proton 
recoil energy of E p ~ 32 MeV. In order to avoid the effects 
of saturation, particularly in the region E p /E n > 1, an 
upper limit of E n — 18 MeV/N has been imposed. 

At low energies the proton recoil is free of saturation ef- 
fects. However, background events arising from 7 and cos- 
mic rays represent a potential contaminant of the E p /E n 
distribution. These events are randomly distributed in 
time, and thus the relative rate increases at low energy 
(Fig. [3]) since Eod~ 2 . As the energy loss in a module is com- 
pletely uncorrelated with the inferred time-of-flight from 
the reaction at the target, E p /E n is not confined below 
1.4. Even if the rejection rate using pulse-shape analysis is 
close to 100 %, any events that remain could mimic a ^n 
signal. In order to reject these events, we have first veri- 
fied that their relative rate at low energy is, as expected for 
a background, independent of the reaction channel (lower 
panel in Fig. [}]). A lower limit on E n was then imposed 
and raised up to the value (11 MeV/N) for which no events 
remain in the ( 14 Be, 12 Be+n) channel, where no ^n should 
be produced. 

The energy gate, E n = 11-18 MeV/N, is shown as the 
shaded area in Fig. |[ As may also be seen in this fig- 
ure, the neutron energy distributions exhibit two compo- 



* The possibility of operating the DEMON photomultipliers at 
lower voltages is being explored for future dedicated experiments. 
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FIG. 3: Data from the reactions C( A Z,X+n) with 14 Be, n Li and 
15 B beams (solid, dashed and dotted lines, respectively). Left- 
top: neutron energy distribution; the shaded area corresponds 
to the energy range used in the present analysis; the percentage 
of events in the hatched area with E p /E n > 1.4 is shown in the 
lower panel as a function of the particle identification parameter 
defined in Eq. ([!]). Right-top: evolution with energy of the ratio 
of 7 and cosmic-ray events to neutrons detected in DEMON. 



nents: the neutrons from the breakup of the projectile (dis- 
tribution centred close to the beam velocity, ~ 30 MeV/N) 
and low energy neutrons evaporated by the excited target- 
like residue. In the case of 15 B, the neutrons arising from 
breakup arc shifted to higher energies due to the higher en- 
ergy of the beam, and therefore the ratio of background to 
neutrons is still relatively high at 11 MeV/N (dotted lines 
in Fig. [|. A limit of E n = 15-18 MeV/N was thus imposed 
for the 15 B data. 



3 Results 

The detection of neutrons produced in the reaction 
C( 14 Be, 12 Be+n) is displayed in Fig. || (histogram); a chan- 
nel in which A n clusters should be absent. We observe that 
the flat distribution predicted for n-p scattering describes 
the data well, except for a small fraction of events at low 
E p /E n . As noted earlier, these correspond to reactions on 
12 C which always generate smaller light outputs p0| . 

The charged fragments produced in the breakup of the 
beam particles were identified using the energy loss (AEsi) 
and residual energy (-Ecsi) signals derived from the tele- 
scope (Fig. |[). One dimensional spectra representing the 
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FIG. 4: Scatter plot of the energies deposited in the Si-Csl tele- 
scope for the reaction C( 14 Be,X+n). Symbols correspond to the 
7 events in Fig. ^ with E p /E n > 1.4. The horizontal band is 
discussed in the text. 
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FIG. 5: Scatter plot and the projections onto both axes of 
the particle identification parameter PID defined in Eq. (|l|) vs 
Ep/En for the data from the reaction C( 14 Be,X+n). The PID 
projection is displayed for all neutron energies. The dotted lines 
correspond to E p /E n — 1.4 and to the region centred on the 
10 Be peak. 
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particle identification were constructed as |27], : 

PID = (AE Si + a) exp {-(E Csl - b) 2 /2c 2 } (1) 

The PID distribution for each beam (left panels in Figs. |||J|) 
exhibits peaks corresponding to isotopes of H, He, Li, Be 
and B. The parameters a, 6, c of Eq. (Q) were adjusted using 
the Be isotopes [|7|, in which the peaks corresponding to 
10,12 Be are well resolved (Figs. The cross-sections for 
the production of 12 ' 11 ' 10 Be from the breakup of 14 Be were 
460 ± 40, 85 ± 15 and 145 ± 20 mb, respectively f27|. 

The E p /E n distributions (upper panels in FigsT^,^) ex- 
hibit a general trend below 1.4: a plateau up to 1 followed 
by a sharp decline, which may be fitted to an exponential 
distribution (dashed line). In the region where A n events 
may be expected to appear, E p /E n > 1.4, some 7 events 
ranging from 1.4 to 2.2 are observed for 14 Be. However, 
in the case of 11 Li, despite the greater number of neutrons 
detected (factor of 2.4) only 4 events appear and lie just 
above threshold, whilst there are no events in the case of 
15 B. Turning to the coincidences with the charged frag- 
ments produced in the breakup, the 7 events produced by 
the 14 Be beam fall within a region centred on 10 Be. In the 
case of the 4 events produced in the reactions with 11 Li no 
similar correlation appears to exist. 

The left panel in Fig. [?] displays in more detail the region 
of the particle identification spectrum for the breakup of 
14 Be into lighter Be isotopes, together with the 7 events 
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FIG. 7: Left: detail of the particle identification spectrum 
around 10 ' 12 Be for the data from the reaction C( 14 Be,X+n); 
the 7 events with E p /E n > 1.4 are denoted by the symbols. 
Right: results of the simulation described in the text for the 
reaction ( 14 Be, 9_12 Be); the shaded histogram is the sum of all 
contributions (lines). 
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FIG. 8: Left: out-of-plane vs in-plane disposition of the DE- 
MON modules as viewed from the target and the location of 
the 6 10 Be+ 3,4 n candidate events (symbols). Right: angular 
distribution. 



for which E p /E n > 1.4. Clearly the resolution in PID does 
not allow us to unambiguously associate all of the observed 
events with a 10 Be fragment. However, the much higher 
cross-section for this channel compared to 9,11 Be does sug- 
gest that this may be the case. 

Reactions in the Si-Csl telescope, which was centred at 
zero degrees ^3|, also affect the PID. If 14 Be traverses the 
target and breaks up in the telescope, the AE$i will corre- 
spond to 14 Be but the Ec s i may take on a range of values 
depending on the interaction point. Such events correspond 
to the horizontal band centred on AE^ ~ 2.5 in Fig. ^. The 
effects of breakup into lighter Be isotopes in the different 
elements of the telescope have been simulated. The results 
(right panel in Fig. |t]) reproduce well the characteristics of 
the observed PID spectrum. Each peak exhibits tails: to 
the left owing to the energy response of the Csl and to the 
right owing to the energy loss of 14 Be in the Csl before 
breakup. Thus, whilst the PID are further broadened by 
the effects of the reactions in the telescope, the 7 events 
in question are most probably associated with a 10 Be frag- 
ment. 

As a first step towards investigating the nature of these 
events we have verified that each corresponds to a well de- 
fined event in both the charged particle and neutron detec- 
tors. As described above, in terms of the charged fragment 
identification (Figs. ^J^), the 7 events lie within the region 
corresponding to 10 Be. One of the events, however, exhibits 
a very low i?csi j such that it is located well away from the 
region of the plot where most of the yield is concentrated. 
This event, with E p /E D = 1.41, is denoted by the open 
symbol in Figs. |],fr| In addition, it is not kinematically con- 
sistent with the breakup of 14 Be into 10 Be+ 3,4 n, as a low 
energy multineutron (11-18 MeV/N) should be associated 
with a high energy 10 Be. In the following the discussion is 
therefore limited to the 6 remaining events. 

Regarding the detection using DEMON, we have checked 
that the 6 candidate events do not cluster in a single de- 



FIG. 9: Pulse-shape discrimination in DEMON for data from 
the reaction C( 14 Be,X+n). Left: integrated slow component 
(Qsiow) vs the total integrated light output (Qtotal) f° r a module; 
the dashed line separates the neutrons (upper group) and 7 or 
cosmic rays (lower group). Right: distance in the previous plot 
to the dashed line as a function of the neutron energy for all 
modules; the symbols correspond to the 6 10 Be+ 3 ' 4 n candidate 
events. 



tector or group of detectors (Fig. g). The operation of each 
of the detectors was then examined in detail — in partic- 
ular the pulse-shape analysis. The standard procedure to 
select the events arising from the neutrons is to compare 
the integrated slow component (Qsiow) of the light output 
with the total integrated output (Qtotal) As displayed 
in Fig. ^, excepting at very low Qtotal, the neutrons form a 
distinct group from the 7 and cosmic-ray events. In order 
to present and compare the 90 modules the data have been 
standardised in terms of the parameter AQ, defined as the 
distance an event lies away from the line separating the two 
groups. The data from the 90 detectors are plotted in Fig. || 
and the 6 events in question are highlighted. As is clearly 
evident all 6 events fall within the region corresponding to 
the detection of neutrons. In order to be sure that any pos- 
sible contamination by 7 and cosmic-ray events (section 2) 
was minimal, the gate shown by the dotted lines in Fig. ffl 
was applied to the data. 

4 Discussion 

The 6 events in question thus exhibit characteristics con- 
sistent with detection of a multineutron cluster from the 
breakup of 14 Be. Any potential sources of events with 
E p /E n > 1.4 not involving the formation of a multineutron 
are now examined in detail. Each scenario must account for 
the fact that the 6 events appear to be produced in asso- 
ciation with 10 Be fragments (17 % of the total yield) while 
no events appear in association with other fragments with 
comparable (H-He-Li, 19 %) or higher ( 12 Be, 47 %) yields. 

We first address the most probable source of events, neu- 
tron pile-up, i.e. the detection for the same event of more 
than one neutron in the same module. Alternative sources, 
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FIG. 10: Data from the reactions C( 14 Be, 12 Be+n), upper, and C( 14 Be, 10 Be+n), lower panels. From left to right: distributions of 
neutron angle, energy and multiplicity. The open symbols are the data, the solid lines are results of the simulations described in 
the text. The histogram in the right-most panel is for a simulation in which 10 7 events were generated; the dashed line is the result 
of a fit to an exponential for E p /E n > 1. 



TABLE 1: Parameters for the simulations shown in Fig. |l0p average neutron multiplicity and momentum width for the projectile- 
like component; average multiplicity, energy slope and angular width for the target-like component; and the pile-up probability 
deduced from the events within the shaded area in Fig. Uu. 



channel 


(M(Be)> 


T [MeV/c] 


(M(C)) 


E [MeV] 6 1/2 [°] 


P 2 „ 


( 14 Be, 12 Be) 


1.4 


85 


1.6 


12 35 


(6± 1) x 10~ 4 


( 14 Be, 10 Be) 


2.8 


140 


4.0 


8 40 


(4± 1) x 10~ 4 



including pile-up of a neutron and a 7-ray, the detection in 
DEMON of light charged particles, and the random back- 
ground which was in principle rejected by the condition 
E n > 11 McV/N, are also reviewed. 

4.1 Neutron pile-up 

Let us consider the possibility of pile-up, or sum events 
|H| , due to the detection of more than one neutron in the 
same module. As the intrinsic neutron detection efficiency 
is moderately low and the detector array highly granular, 
the pile-up of just two neutrons will be the leading con- 
tribution. We define the probability of pile-up leading to 
E p /E n > 1.4 for a given channel X as: 

P 2n (X) = N 2a (X) I N(X) 

where N is the total number of events observed for the 
channel X+n and A^n the number with E p /E n > 1.4. A 



simple estimate of the number of events expected from pile- 
up may be derived from channels in which no events with 
E p /E n > 1.4 arc observed. For example: 

P 2n ( 12 Be) < l/iV( 12 Be) 
P 2n (H-Hc-Li) < 1 /AT(H-He-Li) 

We assume that 2 neutrons are emitted in the 12 Be channel 
and 4 in the 10 Be one. In the latter case, therefore, there are 
3 neutrons available to pile up. For the H-He-Li channel we 
assume emission of at least 4 neutrons. It is also assumed 
that the pile-up distribution in E p /E n is flat and extends 
up to ~ 3, which means that half of these events lie above 
E p /E n = 1.4 (it will be shown below that this considerably 
overestimates the contribution). Therefore: 

A^ 2n ( 10 Be) < N{ w Be) x 3/2P 2n ( 12 Be) = 0.5 
iV 2n ( 10 Be) < N{ w Be) x 1/2 P 2n (H-He-Li) = 0.4 
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At most 0.5 pile-up events are expected, corresponding to 
a probability P 2 n( 10 Be) - 9 x 1CT 4 . 

In this simple estimate we do not consider neutrons from 
the target nucleus, which may differ in energy, angular and 
multiplicity distributions for the different reaction chan- 
nels, or the non-trivial effect of the conditions applied in 
the analysis, such as the limits placed on the energy E n . In 
order to take into account these effects and how they in- 
fluence the pile-up probability in each channel, a complete 
Monte-Carlo simulation has been performed employing the 
code MENATE ||| |§. In the 10 < 12 Be channels, neutrons 
from 1 Be were described by an average multiplicity, with 
a distribution ranging from to 4,2 (respectively), and a 
Lorentzian momentum distribution in the projectile frame 
with the measured width (T) |^3|, |27|. For the neutrons 
from the target nucleus, the simulation also used an aver- 
age multiplicity, with a distribution ranging from up to 6, 
an energy distribution of the form e~ E,1 ^ E ° , and a Gaussian 
angular distribution with a half- width #1/2- 

The different parameters were chosen so as to reproduce, 
in each channel, the measured neutron energy, angular and 
multiplicity distributions. These distributions (open sym- 
bols in Fig. [To]) display the main differences between the 
neutrons emitted in the two channels: (i) the target nucleus 
contribution for 10 Be is larger and exhibits a smaller slope; 
(ii) neutrons in coincidence with 10 Be exhibit a broader 
momentum distribution; (iii) the multiplicity is higher for 
10 Be. The simulation (solid lines), using the parameters 
listed in Table [|, provides a good overall description of all 
these features. 

The histograms in Fig. |l^ represent the E p /E n distribu- 
tions obtained after generating 10 7 events and applying the 
same analysis conditions as those applied to the data from 
the experiment. It should also be noted that the simulation 
includes the pile-up arising from more than two neutrons 
interacting in a module. The events resulting from pile-up 
are clearly evident for E p /E n above 1.4. Furthermore, as 
expected for pile-up, the number of events decreases with 
increasing E p /E n . Interestingly, P 2n for the 10 Be and 12 Be 
channels are comparable (Table [t]) despite the higher mul- 
tiplicity for the former — this is most probably due to the 
enhanced forward focussing of the halo neutrons in the lat- 
ter channel. The estimates lead to iV2 n ( 12 Be) = 0.8, in 
agreement with the experiment, and 7V2 n ( 10 Be) = 0.2. As- 
suming P 2n ~ 5 x 10~ 4 as an average value for the other 
two beams, iV2 n ~ 3.3 for n Li and A^n ~ 0.3 for 15 B, also 
in agreement with experiment. 

Another estimate of pile-up, which does not rely on sim- 
ulations, may be made. By selecting multiplicity 2 events, 
the number of events in which two neutrons interact in the 
same detector can be deduced from the relative-angle dis- 
tribution. Importantly, this approach includes explicitly 
all the effects arising from the different conditions applied, 
as the experimental distributions (Fig. |ll"| ) have been ex- 
tracted under analysis conditions equivalent to those used 
to derive the results in Figs. The distributions have 
been fitted according to a Gaussian lineshape for 9 nn > 2° 




10 20 10 20 10 20 



FIG. 11: Distribution of relative angle for neutron pairs from 
the reactions C( A Z,X+2n) with 14 Be and n Li beams, selected 
with conditions equivalent to those used in Figs. |^^. The solid 
line is the result of a fit to a Gaussian, which is extrapolated to 
6> nn < 2° (dotted line). 



TABLE 2: Comparison of the number of events with E p /E n > 
1.4 for each channel (1 st column) with the estimated number 
of events expected from pile-up. The methods are based on the 
absence of events in the 12 Be channel (2 nd column), Monte-Carlo 
simulations (3 rd column), and the relative-angle distribution of 
n-n pairs (4 th column). 



channel 


Jv 2n 


iv 2n 


Ar (sim) 
Jv 2n 


Ar(nn) 
iv 2n 


( U Li,X) 


4 


<6.0 


-3.3 


<7.0 


( 15 B,X) 





<0.5 


-0.3 


<0.9 


( 14 Be, 12 Be) 







0.8 


<1.2 


( 14 Bc, 10 Be) 


6 


<0.5 


0.2 


<0.8 



(solid lines) . The number of pairs that may be expected to 
be detected in a single module can now be estimated from 
the extrapolated fit (a single module on average subtended 
1°). Using the value of the fit at 1° and assuming that 
half the pairs result in E p /E n > 1.4, the estimated pile-up 
probability is: 

P2n( 10 Be) - 12 x 10 -4 

It should be stressed that this estimate represents a 
very conservative upper limit. Even if the uncertainty of 
the extrapolated fit at 1° is of the order of 20 % for the 
10 ' 12 Be channels (Fig. [ll]), the number of pile-up events 
with E p /E n > 1.4, as indicated by the simulations, is much 
less (1/3-1/6) than half the number of double hits. This 
method thus suggests that iV 2n ( 10 Be) < 0.8, iV 2n ( 12 Be) < 
1.2, N 2n < 7.0 for n Li and N 2n < 0.9 for 15 B. These results 
are in line, except in the case of 10 Be, with the number of 
observed events. A summary of the different estimates of 
the number of pile-up events is given in Table |[ 
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4.2 Other sources of background 

A number of other sources of background may be postu- 
lated to account for the excess of events with E p /E n > 1.4. 
Above all, any such processes must account for the fact 
that these events appear to be correlated with a 10 Be frag- 
ment. One possible process that can be eliminated is the 
production of a large number of 7-rays associated with the 
( 14 Be, 10 Be) reaction channel. Pile-up of a neutron with a 
7-ray could conceivably produce events with E p /E n > 1.4, 
but the flight-time for such events would be that of the 
7-ray and therefore the event would be rejected. 

Other processes involving the detection of particles not 
arising from reactions of the 1 Be beam particle cannot be 
at the origin of the events in question, as they would not 
give rise to a correlation with a particular channel. Light 
charged particles, for example, if present in the beam and 
yet not triggering the Si beam-identification detector or be- 
ing detected in the Si-Csl telescope, should show no corre- 
lation with any particular charged fragment. In addition, 
such events would be concentrated very close to the beam 
axis, which is not the case for the 6 events observed here 
(Fig. |). 

Finally, for background events in DEMON leading to 
E p /E n > 1.4 for an apparent E n < 11 MeV/N (section 2), 
we have considered the correlation with the charged frag- 
ment identification (PID). The result is shown in the lower 
panel of Fig. || for the 14 Be beam data. The rate of events, 
~ 5 %, is independent of the PID, suggesting that such 
a contamination cannot be at the origin of the observed 
events. 



4.3 Kinematical constraints and properties 

We have analysed the kinematics of the 6 events that 
are consistent with the detection of a bound 3,4 n cluster. 
Breakup reactions in the 10 Be channel leading to no multi- 
neutrons, a trineutron, or a tetraneutron are: 



10 Be+4n 

10 Bc+ 3 n+n 

10 Be+ 4 n 



(2) 
(3) 
(4) 



with five, three and two particles in the final state, respec- 
tively. In the average-beam-velocity frame, therefore, the 
relative angle (8*) between the neutron cluster detected and 
the 10 Be fragment should increase with decreasing particle 
number, up to a value close to 180° for reaction (||). 

The experimental distributions of relative angle are 
shown in Fig. |l2|: the solid line corresponds to reaction , 
while the dashed line is for the 12 Be channel, which is sim- 
ilar to reaction (|^) in that it has three particles in the final 
state. Indeed, both distributions are centred at backward 
angles. This is due to the selection of neutrons with en- 
ergies lower than the beam, as in the beam-velocity frame 
these always correspond to neutrons emitted backwards. 
However, the expected trend is observed when going from 




[deg] 



180 20 40 20 40 

E [MeV/N] m[MeV] 



FIG. 12: Distributions for the 6 multineutron candidate events 
(histograms). Left: relative angle between the Be fragment 
and the neutron cluster in the average-beam-velocity frame, 
compared to the data for 10 Be+n (solid line) and to those for 
12 Be+n (dashed line); both distributions have been normalised 
to 6 events. Middle-right: total kinetic energy and invariant 
mass with respect to the threshold for the 10 Be+ 4 n system; the 
dashed lines are the result of a simulation in which the input 
distribution (dotted line) has been filtered by the analysis con- 
ditions used in Fig. pi 



a five- to a three-particle final state, as demonstrated by 
the increase in the average value from (8*) — 139° to 146°. 

The relative angle between 10 Be and the neutron cluster 
for the 6 events is biased more closely to 180°, as expected 
for a back-to-back decay. The average value is (8*) = 157°. 
This argues in favour of reaction (|J), which suggests that 
the most probable scenario for the events observed is the 
production of a tetraneutron in the breakup of 14 Be. 

We have subsequently considered the analysis of the en- 
ergies assuming that 10 Be+ 4 n breakup took place. As the 
tetraneutron energy has been restricted to 11-18 MeV/N 
and the 14 Be beam energy is around 30 MeV/N, 10 Be frag- 
ments in coincidence with 4 n should have energies per nu- 
cleon higher than the beam. The resulting total kinetic en- 
ergy (Fig. middle panel) corresponds well to the beam 
energy. Moreover, the reconstructed 14 Be* invariant mass 
might show strength at energies just above the 10 Be+ 4 n 
threshold (S^n = 5 MeV), as was observed for 12 Be+2n 
p3fl . The invariant mass distribution (Fig. [l2|, right panel), 
however, is located far above the threshold. This is related 
to the low-energy gating of the neutral fragment since, as 
noted above, the corresponding 10 Be fragments have high 
energies and thus the relative energy of the two fragments 
is high. This interpretation is confirmed by the dashed lines 
in Fig. [12], middle-right panels, which correspond to the re- 
sults of a simulation of reaction (jij) with a 14 Be invariant 
mass with strength at low energy (dotted line) and the low- 
energy gate on 4 n. They are in good agreement with the 
experimental distributions. 

The events observed here may be used to provide esti- 
mates of the momentum content, lifetime and production 
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cross-section of the postulated tetraneutron. A fit of the 
angular distribution using a Lorentzian lineshape (Fig. |^) 
leads to a momentum width of V = 90 ± 60 MeV/c. At 
best, it can be concluded that this is compatible with the 
widths measured for the single neutrons in the breakup of 
14 Be §§. 

Concerning the lifetime, the events observed suggest that 
4 n is particle stable. However, it will be unstable against 
f3 decay because the binding energy is limited to 3.1 MeV 
(section 1), which implies p3| : 

m( 4 n) - {m( 3 H) + m n } > 6.2 MeV 
m( 4 n) - m( 4 H) > 3.2 MeV 

The average flight time of the 6 events from the target to 
DEMON is ~ 100 ns. This indicates that the lifetime must 
be of this order or longer. 

The conditions applied in the analysis presented here 
make an estimate of the total 4 n production cross-section 
rather difficult. Nonetheless, if we assume that the various 
gates affect the number of neutrons and tetraneutrons in a 
similar manner, we can scale the cross-section measured for 
the production of 10 Be [^7) by the relative yield observed 
and obtain er( 4 n) ~ 1 mb. 

5 Summary 

A new approach to the production and detection of multi- 
neutron clusters has been presented. The technique is based 
on the breakup of energetic beams of very neutron-rich nu- 
clei and the subsequent detection of the multineutron clus- 
ter in liquid scintillator modules. The detection in the scin- 
tillator is accomplished via the measurement of the energy 
of the recoiling proton (E p ). This is then compared with 
the energy derived from the flight time (E a ), possible multi- 
neutron events being associated with values of E p > E n . 

The method has been applied to data from the breakup 
of 11 Li, 14 Be and 15 B by a C target. In the case of the 
14 Be beam, some 6 events have been observed with charac- 
teristics consistent with the production and detection of a 
multineutron cluster. Special care has been taken to esti- 
mate the effects of pile-up. Three independent approaches 
were applied and it was concluded that at most pile-up may 
account for some 10 % of the observed signal. As discussed, 
the most probable scenario was concluded to be the forma- 
tion of a bound tetraneutron in coincidence with 10 Bc. 

The confirmation of the events observed here with a 
higher intensity 14 Be beam and an improved charged parti- 
cle identification system, and the search for similar events 
in the breakup of 8 He, appear to be the most straightfor- 
ward steps to take in the near future. The saturation effects 
encountered with DEMON at high light outputs should be 
reduced considerably by lowering the beam energy and the 
high voltage applied to the photomultipliers. In the longer 
term, searches for heavier multineutron clusters may be 
conducted when more neutron-rich beams become available 
at intensities beyond 10 2 pps. 



Theoretically, multineutron clusters may provide a strin- 
gent test of our basic understanding of the nuclear inter- 
action. If the observed events are confirmed, the challenge 
will become one of describing how neutron clusters may be 
bound. Possible avenues to such an end include advances in 
few-body theory, changes in parameters describing the N-N 
potential, or the introduction of new terms or many-body 
forces that occur in proton-free environments. 
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